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Abstract 

We analyze the new equation of motion for the damped osciUator. It differs from the standard 
one by a damping term which is nonlocal in time and hence it gives rise to a system with memory. 
Both classical and quantum analysis is performed. The characteristic feature of this nonlocal 
system is that it breaks local composition low for the classical Hamiltonian dynamics and the 
corresponding quantum propagator. 
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1 Introduction 

The damped harmonic osciUator is one of the simplest quantum systems displaying the dissipation 
of energy. Moreover, it is of great physical importance and has found many apphcations especially 
in quantum optics. For example it plays a central role in the quantum theory of lasers and masers 
[I] [51 [3] . The standard Hilbert space formulation of Quantum Mechanics gives rise to the unitary 
evolution with no room for dissipative behavior. The usual approach to include dissipation is the 
quantum theory of open systems [11 [6l [71 [8] . In this approach the dynamics of a quantum system is 
no longer unitary but it is defined by a semigroup of completely positive maps in the space of density 
operators 0. 

Another strategy to describe dissipative quantum systems is based on the old idea of Bateman [TDj . 
Bateman has shown that in order to apply the standard canonical formalism of classical mechanics 
to dissipative and non-Hamiltonian systems, one can double the number of degrees of freedom, so 
as to deal with an effective isolated classical Hamiltonian system. In this approach the new degrees 
of freedom may be interpreted as variables describing a reservoir. Applying this idea to a damped 
harmonic oscillator one obtains a pair of damped oscillators (so called Bateman's dual system): a 
primary one and its time reversed image. The Bateman dual Hamiltonian has been rediscovered by 
Morse and Feshbach [11] and Bopp [12] and the detailed quantum mechanical analysis was performed 
by Feshbach and Tikochinski [13j . The quantum Bateman system was then analyzed by many authors 
(see for the historical review [14j). Recently it was investigated in connection with quantum field 
theory and quantum groups [15] (see also [El [17]). Finally, the detailed analysis of spectral properties 
of the Bateman system was recently performed in [TB] [TH] (see also [10] )• For other approaches to 
quantization of damped oscillator see also recent papers [lH [22l [23] . 
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In the present paper we propose a different approach. Instead of analyzing the standard equation 
of motion for the damped osciUator 

x + + = , (1) 

where 7 is a positive damping constant, and a; is a frequency of the undamped motion, we shall study 
a modified equation 

i + 27tanh(7(t-io))a;+t^^a; = , (2) 

where <o is an arbitrary constant. It is clear that in the limit ^^t oo\ formula ^ reproduces 
the standard equation ([T|). Therefore, both dynamics may differ for short times only, that is, for 't' 
satisfying 

t <l/7 • 

What is the origin of ([2])? It turns out that both ([1]) and ([2]) give the same frequency of oscillations 
of the damped motior0 

n = - 7' ■ (3) 

Moreover, there is no other linear equation of motion with this property. Interestingly, equation ((2]) 
introduces additional parameter Iq. It turns out that to is responsible for an effective memory, that 
is, in addition to a 'local time evolution governed by ^ 'remembers' about its history at to. It 
means that the initial value problem for ([2]) has to be supplemented by initial data x{to) = xo and 
x{to) = vq. There is a crucial difference between ([Ij and ([2]): the initial data for the standard 'local' 
equation may be provided for arbitrary initial moment of time. It is no longer true for the modified 
'nonlocal' equation: now the initial moment necessarily equals ^o- Changing initial moment to one 
changes the corresponding equation ^ and hence the whole evolution of our damped system. This 
property is responsible for a peculiar feature of ([2]). The composition law for classical dynamics is no 
longer true, that is, evolution from to to ti, and then from ti to t2 differs from the direct evolution 
from to to t2- 

Now, the quantum system corresponding to 12]) may be easily constructed. Our damped system 
turns out to be governed by nonlocal Hamiltonian and hence the standard composition law for quantum 
propagators 

K{x2,t2;xo,to) = J dxiK{x2,t2; xi,ti)K{xi,ti]Xo, to) , (4) 

does not hold. We stress that the local composition law (jl]) is satisfied by the quantum propagator 
derived from the standard equation ([1]) via well known Caldirola-Kanai approach ^24, ,25^ . The vio- 
lation of (jll defines the basic difference between two descriptions based on ^ and ([2]), respectively. 
We hope that our analysis sheds new light on this old problem. 

2 Hamiltonian formulation 

There were many attempts for hamiltonian formulation of Newton equation for systems with one 
degree of freedom (see e.g. [351 [221 123 [Ml US] ) ■ K is well known that ([1]) may be derived from the so 
called Caldirola-Kanai Hamiltonian [24l [25] 

Hix,p,t) = ^moio'x'e'-" + -^p'e-'''' . (5) 
2 zmo 

-"^ Throughout the paper we assume i^; > 7. 
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It is clear that Caldirola-Kanai Hamiltonian may be considered as the standard Hamiltonian of the 
harmonic oscillator with time-dependent mass given by 

m{t) = TOoe^'^* . (6) 

Consider now a family of Hamiltonians with an arbitrary time-dependent mass m(t): 

= ^ + (7) 
The corresponding Newton equation reads as follows 

X+^^^^X + LJ^X = . (8) 

m[tj 

Let us assume that in the asymptotic regime ^jt — )■ oo', equation ([8|) reproduces the standard equation 
(d]), that is 

2 m[t) 

To solve (|5]) one represents complexified x{t) as follows 

x{t) = X{t) e"^^*^ , (10) 

and obtains the following equations for time-dependent real quantities A, (f> and k: 

X- X<P^ +2kX + uj^X = 0, (11) 
(jiA + 20A + 2k</)A = 0, (12) 

where we omitted explicit time-dependence. Now, we make our basic assumption that = 17 with 
defined in ([3]). Equation p2)) reduces then to 

A + kA = , 

and hence equation (fTT|) implies 

o2 2 2 

ll = UJ — K — K , 

which together with ^ gives finally 

• I 2 2 

K + K =7 . 
The above equation has the following solutions 

k(<) = 7tanh(7(i - to)) or K{t) = ±7, (13) 

with to being an integration constant. Note, however, that due to ^ the solution K{t) = —7 is not 
acceptable. It is clear that K{t) = 7 corresponds to the Caldirola-Kanai model with m(t) = e^'''*, 
whereas K(t) — 7tanh(7(i — to)) gives rise to our new equation ^ described by a Hamiltonian with 
time-dependent mass 

m(t) = mo cosh^(7(t - to)) , (14) 

satisfying m(to) — niQ. We may therefore interpret (l2|) as a Newton equation for a particle with 
constant mass mo under the action of the time dependent nonlocal force 

F{t; to) = — 2mo7 tanh(7(t — to))i — nioui'^x , 

or as a Newton equation for a particle with time-dependent nonlocal mass (I14p under the action of 
perfectly local force ^niQUj^x. 
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3 Classical dynamics 

Consider now classical Hamiltonian dynamics generated by ([8]) with m(t) defined in (jl4p . It is given 

by 



where the 2x2 matrix 



reads as follows 



A new 

XX 



x{t) 
Pit) 



cos{n{t — to)) 

cosh(7(t - to)) 



A"™(t;to) 



Pq 



A new A new 
XX xp 

A new A new 
px pp 



A new 

xp 



1 sm{n{t - to)) 
moil cosh(7(t — to)) 



and 



A new 

px 



A new 
■^^pp 



'Too 



n cosli(7(i -to)) sin{fl{t - to)) + 7 sinli(7(t -to)) cos(ri(t - to)) 



7 . 



Q cosh(7(t — to)) cos(51(t — to)) — — sinh(7(t — to)) sin(il(t — to)) 



The crucial feature of (|T5l) is that it breaks the local composition law 

A(t2;to) = A(t2;ti)oA(ti;to) , 



(15) 



(16) 

(17) 
(18) 

(19) 

for to < ti < t2, that is, formula (HH) does not hold for A = A"™. 

Let us compare our new evolution (jlSp with the classical dynamics of the Caldirola-Kanai system. 
Note, that we may trivially change the formula ^ for time-dependent mass 



TO(t) = TOoe^'^f*"*") , 



(20) 



without affecting the equation of motion. Now, these two models contain additional parameter to. 
One easily finds for the corresponding Caldirola-Kanai dynamics A*-^^(t;to): 



^CK 

XX 



CK 

xp 



cos(fi(t - to)) + ^ sin(f7(t - to)) 
^sin(f^(t-to))e-^(*'*°), 

TOoSi 



-7(t-*o) 



CK 



a: 



^CK 

PP 



-mon{ 1 + ^) sm{n{t - to))e^(*-*°\ 



cos(fi(t - to)) - ^ sin(f7(t - to)) 



o7(t-to) 



One may check that A*-^^ satisfies 

ACK(t2;to) = A^'^it2;ti)oA''^{tv,to) , 



(21) 
(22) 

(23) 
(24) 

(25) 



for to < ti < t2. Note, that computing A*-^^(t2,ti) one applies (|2ip - ([24|) with to replaced by ti, and 
Too replaced by nii := m(ti) = moe^'*'^*i~*°-'. Then formula easily follows as a simple consequence 
of e°+^ = e°e''. It is no longer true if we replace (PU]) by (fH|). 
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Finally, let us compare both approaches in the simplified situation corresponding to pure dissipa- 
tion, that is, uj = 0. In this case the corresponding formula for A"°'"(t; to) considerably simplifies 

A^it-M = I (™o7)-^tanh(7(i-io)) ^ (26) 

However, Aq°* still breaks p^ . In the case of Caldirola-Kanai dynamics one finds 
A--ifM ^ ( I (™o7)-sinh(,(i-io))e--(-«) ^ _ 

Note, that in a asymptotic regime — *■ 00 one obtains 

Anew ^ f 1 {mair^ \ aCK ^ f 1 (2mo7)"^ \ ^r,oN 

■^*-asympt 1 1 / ' asympt \ 1 / ' \ / 

that is, both dynamics have different asymptotic states: (xq +Po/(™o7):Po) smd [xq +Po/(2too7)jPo) 
for the Caldirola-Kanai system. 

4 Quantum propagator 

Consider now a quantum propagator corresponding to Hamiltonian ([7]) with m{t) given by (|14p . Since 
the Hamiltonian is quadratic in one easily finds for the propagator (see e.g. |32p 

K{x, t; xo, to) = N{t - to) A{x, t; xqM) ■ B{x, t; to) , (29) 

where 



A{x,t]Xo,to) = exp . /^f'f^ ^ (xq + x'^) cos(n(t — to)) ^ 2a;oxcosh(7(< — io)) 

n cos{n{t - to)) sinh(7(i - to)) 



2hfi\n{n{t-to)) 

4- + \ J * 2 sinh(7(t ~ to)) 

B[x,t]to) — exp — mox 



2h sin(f^(t - to)) 

7sin(r!(t-to))cosh(7(t-to))]} , (30) 



and the normalization factor is given by 



y 27ri/ism(ii(t — to)) 

Note, that when 7 = the factor B equals to 1 and one recovers standard propagator for the harmonic 
oscillator [5D]. On the other hand when a; = 0, that is, in a pure damping case, formula (^5)1 reduces 
to 



^ 

.2fttanh(7(t - to)) / 

with 



/r„(...;.„.M^A'.(i-i.)exp{|I J° !;-;°' , } . (32) 
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For 7 = formula ([5^ reproduces quantum propagator for a free particle [50]. Interestingly, the 
formula ([32]) for a quantum propagator for a particle undergoing linear damping was recently derived 
by Kochan [3T] by using purely geometrical considerations. Here we showed that it may be derived 
via standard methods by taking suitable nonlocal Hamiltonian with memory. In [31j the formula for a 
full propagator of the damped oscillator was also presented which however differs from (P^ . Actually, 
Kochan propagator is given by 

^Kochan(a;, t; xo,to) = N{t - to) A{x, t; xq, io) , (34) 

and therefore differs from (^5)1 only by a factor B{x,t;to). Both propagators K and /^Kochan produce 
the same limits for 7 = (harmonic oscillator) and for (7 = & = 0) (free particle). Observe that 
B{x, t; to) produces only an additional phase factor at the final position 'x\ Therefore both formulae 
give the same results for a mean position of a particle: starting from the initial wave function ^{xq, tg) 
one finds 

00 

ip{x,t;to)^ J K{x,t;xo,to)i{jixo,to)dxo , (35) 
—00 

and 

00 

V'Kochan(a;,t;<o) = / KKochiin{x,t;xo,to)i>ixo,tQ)dxo = — — — — V^(x,t;io) , (36) 
J B(x,t;to) 

— 00 

and hence 

00 00 
{x{t-to)) = J x\ij{x,t;to)\'^dx ^ j x\iiKoc\vs.n{x,t;to)\^dx . (37) 

— 00 — CX3 

It is easy to show that {x{t;to)) satisfies our modified nonlocal equation ([2]): 

d 

-j^(x{t;to)) + 2-ftanh{-f{t - to))-{x{t;to)) + Lo^x{t;to)) ^ . (38) 

In the Appendix we show that this additional phase factor B{x, t; to) may be eliminated by a suitable 
canonical transformation. The price we pay for a simpler formula for a propagator is the much more 
complex form of the Hamiltonian which apart from the standard quadratic terms x^ and contains 
also a mixed term xp (see Appendix). 

Let us compare the above formulae with the propagator for the Caldirola-Kanai system. One 
easily finds the CK propagator corresponding to ■m{t) given by ((20)) 

KcK{x,t;xo,to) = iVcK(i-io)exp{^mo[7(a;^e"^*° -x2g27t^ 

+ ncomt-to)){.y^-^.^e^^^) - \,,J_,.f ^oe^^''''"^]}, (39) 



where the normalization factor 



NcKit to) = J ' ' , ■ (40) 
y 27rin,sm(i2(t — to)) 
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Again ([5^ correctly reproduces both the propagator for the harmonic oscillator (for 7 = 0) and the 
propagator for the free particle (7 = and u = 0). However, for pure damping both propagators are 
different. 

There is a crucial difference between K and Kqk- The Caldirola-Kanai propagator satisfies local 
composition law 



(41) 



ft is no longer true for K defined in ([^5]) . Therefore, as in the classical case memory term encoded 
into dynamical equation is responsible for the breaking of the local composition law for the quantum 
evolution. 

Finally, let us consider a damped (without oscillations) evolution of the initial Gaussian wave 
packet 

Mx; to) = Ce-(^-^»)'/'^+*P«^/'' . (42) 
One easily finds for the probability density p{x,t;to) = \4>{x,t;tQ)\'^: 



p{x,t;to) 



\C\' 



: exp 



(x-x(t)f 



v/cr(t - to) /a 
where x{t) is given by the classical formula 

tanh(7(t — to)) 



cr{t - to) 



x(t) = xo +po- 



mo7 



and 



^(t -to) = a- 



1 



/itanh(7(t — to)) 
mou^ 



(43) 



(44) 



(45) 



It shows that the presence of the damping term modifies the way the initial way function spreads in 
time. Note, that in the asymptotic regime one obtains 



Pasympt ('^' ; ^o) 



\/ Casympt/f 

with the classical formula for the asymptotic position 



exp 



-2 



{x ^asympt) 
"^asympt 



f-asympt ^ ^^0 



mo7 



and asymptotic dispersion 



^asyinpt — ^ 



1 



(46) 



(47) 



(48) 



^TOo7cr/ 

It is therefore clear that the damping prevents the wave function from the total spreading. An 
initial Gaussian probability density poix^to) cannot 'spread more' than Pasympt(a;; to)- A similar 
phenomenon appears for Caldirola-Kanai dynamics with different asymptotic state 



Xo 



Po 



-asympt ■ ^^^^ 

but with the same asymptotic dispersion (Jasympt given by 



(49) 
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5 Conclusions 



We performed an analysis of the new equation of motion for the damped osciUator ([2]) ■ It differs from 
the standard one by a damping term - 2"ftainh('j{t — to))x - which is non-local in time and nonlinear in 
the damping constant 7. The new parameter introduces effective memory. For long time behavior 
t ^ 1/7 one recovers standard equation with a damping term - 2jx. Both classical and quantum anal- 
ysis is performed. The characteristic feature of this nonlocal system is that it breaks local composition 
law for the classical Hamiltonian dynamics and the corresponding quantum propagator. Interestingly, 
the same propagator was recently derived in [31| . Without referring to Hamiltonian formulation and 
using purely geometric methods the author of [31 derived the corresponding propagator starting 
from the standard equation ([1]). We have shown that the corresponding classical limit leads to the 
new damping equation Therefore, we conclude that the quantum nonlocal propagator derived in 
[31j does originate in nonlocal equation Finally, it was shown that the purely damped behavior 
modifies well known property of the free quantum dynamics leading to the perfect spreading of the 
initial wave function. It is no longer true when the dissipation is present. Now, instead of the perfect 
spreading there is an asymptotic state giving rise to the asymptotic probability density (|46p . 

Appendix 

Using a general formula for a propagator corresponding to quadratic Hamiltonian (see e.g. [32]) one 
easily finds that 

H{x,p,t) = ^^(t)p^ + :y{t)xp + ^Xit)x\ (A.l) 

with 

,j,{t) = ^ (A.2) 

mo cosh {j{t — to}) 

[27 tanh(7(t~to)) ^ 1 - [1 + tanh^(7ft - tp))] cos^(r!(t - tp)) ] 
^ tan^t-tp)) + .in-mt-to)) / ' ^^''^ 

produces p4|) . that is, it annihilates additional phase factor B{x,t;to) out of ([29]). 
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